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Abstract. We consider the hydrodynamics of an incompressible fluid on a 2D 
periodic domain. There exists a family of stationary solutions with vorticity 
given by Q* = Qcos(p-x)+/3sin(p-x). This situation can be approximated as 
a structure preserving finite dimensional Hamiltonian system by a truncation 
introduced by |24l I26| or by the more standard Galerkin style finite element 
method. We use these two truncations to analyse the linear stability of these 
solutions and analytical and numerical results are compared. Following the 
methods used by El the problem is divided into subsystems and we prove 
that most subsystems are linearly stable. We derive a sufficient condition for 
a subsystem to be linearly unstable and derive an explicit lower bound for the 
associated real eigenvalues independent of the truncation size N. Then we 
show that the corresponding eigenvectors are in . This together with known 
stability results for the 2D periodic Euler equations allows us to conclude 
that most of these stationary solutions are nonlinearly unstable. We confirm 
our results with a numerical computation of the spectrum for a large, finite 
truncation. Finally we discuss the essential spectrum of the full problem as 
the limit of the truncated problem. 


1. Introduction 

In terms of the vorticity ri(x, i) : (T^ x R’*') —>■ R, the 2D incompressible Euler 
equations are (see [I] Appendix 2 for an overview) 


( 1 . 1 ) 


dfl dH. dH. dui du2 


Here x = {xi,X 2 )'^ and ui, U 2 are the velocity components in the xi and X 2 
directions respectively. We impose periodic boundary conditions rt{Tr,X 2 ,t) = 
Q.{—TT,X 2 ,t) and Tl(xi,TT,t) = Tl(xi,—'K,t). There is a family of stationary solu¬ 
tions given by 


(1.2) = a cos(p • x) -I- /3 sin(p • x) 

for a, /3 G R and p G 

The study of stability of certain solutions of the planar Euler equations was 
initiated by the seminal work by [3] on the Lie-Poisson structure of the Euler 
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equations where he invented the Energy-Casimir method to prove stability. This is 
revisited in [2], in particular Section II.4. Arnol’d discusses a slightly more general 
problem where the torus has dimensions X x 2 tt and p = (0, 1)^, and shows that 
this solution is non-linearly stable when X < 2?!. In m it was shown that any 
equilibrium with p = (0,1)^ and X > 27r is linearly unstable for the viscous 
problem using linear stability analysis and infinite continued fractions. That paper 
also shows linear instability for p = (0,m)^, m > 1 and any X. The linear 
instability of the inviscid problem for p = (0,m)^, m > 1 was proved in and 
discussed in [7]. Under the condition ^ for any positive integers 

mi,m 2 , it was shown in |10) that the steady state for p = (0,to)^ is nonlinearly 
unstable. 

In [17] it is shown how to block-diagonalise the linearisation about the equi¬ 
librium with general p into so-called ‘classes’, and using this approach he again 
showed that p = (0,1)^ is Lyapunov stable. This is used in (T^j where the essen¬ 
tial and discrete spectrum of the linearisation of (EB) are studied at the steady 
state EH). They studied the full infinite system, approaching the problem from 
a functional analytic perspective. They found an upper bound on the number of 
non-imaginary isolated eigenvalues, and described the essential spectrum. Further¬ 
more, they showed that the spectrum of the linearised operator is the union of the 
spectrum coming from each of the classes from m and that the spectral mapping 
theorem holds for the Euler operator linearised about 17* (Theorem 2 in [l6|l. 

In the viscous problem solutions * cos{mx) are called bar states in [4]. They 
show that for non-zero viscosity u, and m = 1 these bar states are ‘quasi-stationary’, 
in that they decay on a slow timescale depending on the viscosity. 

We combine the block-diagonalisation used by m with the structure preserving 
finite-dimensional sine truncation [24] and the Galerkin (see [18]) finite element 
truncation to prove that for a large class of p the stationary solutions EB are 
nonlinearly unstable. Zeitlin’s sine truncation leads to a finite dimensional Poisson 
structure and the Hamiltonian structure of the original PDE and its Casimirs are 
preserved in this finite-dimensional truncation. The Galerkin truncation does not 
preserve these Casimirs. See [2] and [T^ for a discussion of the use of Poisson 
brackets in hydrodynamics. The theoretical background of Zeitlin’s sine truncation 
(and a related truncation for a spherical domain) is discussed in [12], and the 
concept of a “limit” of this algebra is discussed in depth in [B] . 

In Sectionj^ the problem and the associated notation are introduced. The system 
is first decomposed into Fourier modes which are described by a non-canonical 
infinite dimensional Hamiltonian system. Then truncation is taken to reduce to a 
finite-mode approximation. We linearise around the steady state, which decouples 
the problem into subsystems. 

In Section |3l we reproduce the “stable disc theorem” from [17] in the truncated 
setting. This theorem states that for classes whose mode numbers a satisfy |a| > |p| 
the spectrum is stable. Thus most class subsystems do not contribute unstable 
modes to the spectrum of the full operator. Then we prove our “unstable disc 
theorem” o which states that if exactly one mode number of a given class is inside 
the unstable disc then for sufficiently large N there is a positive real eigenvalue. In 
our fundamental Theorem 13.51 we show that with certain additional assumptions 
this real eigenvalue is bounded away from zero when N ^ oo. Furthermore in 
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Lemma lS.BI we show that the corresponding eigenvector of the infinite dimensional 
system is in . 

Section 21 provides the main Theorem id. 31 demonstrating non-linear instability of 
the stationary solution m for all choices of p but a few exceptions. In Lemma 14.21 
we establish when the conditions for the lower bound needed in Section [3] are met. 
Zeitlin’s truncation requires some care when proving results for both stable and 
unstable classes. Specifically it is not clear that the intersections between our 
classes and the disc |a| < |p| behave in the way we expect. In Lemma HTT] we show 
that for most choices of p, there is an appropriate truncation size N to control the 
Zeitlin truncation so that Theorem 13.51 can be applied. The other cases of p can be 
treated using the Galerkin trunction. 

Section [5] provides some numerical results. The numerical efficiency and accu¬ 
racy of Zeitlin’s truncation is compared favourably to the Galerkin truncation. A 
connection is made between the nature of the subsystems and the number and type 
of non-imaginary eigenvalues. A discussion of the number of non-imaginary eigen¬ 
values is included, and the accuracy of our calculated lower bound is assessed. A 
brief section on the pure imaginary spectrum of our finite mode systems is included, 
replicating the results in [16] via a very different method. We can show that the 
pure imaginary spectrum of our finite dimensional approximation approaches the 
essential spectrum of the full system. As a result we can naturally define a density 
of eigenvalues in the essential spectrum. 


2. VoRTiciTY Evolution in Fourier Space, Truncation, Linearisation 


2.1. Hamiltonian Formulation. The stream function iL is defined through its 
relation to the fluid velocities by 


( 2 . 1 ) 


Ui 


”*”0—’ ^2 — —• 

0x2 UXi 


The relationship between the stream function and the vorticity is 


( 2 . 2 ) 


n = 


and hence the PDE can be written as 

on _ dn d'^ on aT 

dt dxi dx2 dx2 dxi 

For a fixed p S and T G IR we wish to analyse the steady state 11* = 
a cos(p • x) -I- /3 sin(p • x). 

Note that we can write H* = 2rcos(p • x -f 0), where 9 = ±tan“^ (^) 

\/(y2j^Q2 

r = —2 -■ The signs of 9 and T will depend on the signs of 0 . and /3. If 

0 = 0, then take 6* = |. Define c = j^P, so H* = 2r cos(p • (x -|- c)). Thus 
by taking the translation x = x -|- c we can instead just consider the steady state 
n* = 2r cos(p • x) by a change of origin. Therefore for the remainder of this paper 
we drop the tilde and simply consider the steady state 

D* = 2r cos(p • x). 

Expand D into a Fourier series with coefficients Wk(t) as 

ll(x,t) = ^ Wk(l)e*'"'’" 
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and combine (12.21) and (1^ . Then the Fourier coefficients are governed by the 
ODEs 

. , . / \ k X 1 

(2.4) Wk(t) = 2^ —yi^W-lUk+l 

iez2\{o} ' ' 


(where x x y = Xiy 2 — X 2 yi for x,y G ]R^, and Wk := ^(^k))- The condition 
Wk = wlk is necessary for £7 to be real. 

Define the ‘ideal fluid’ Poisson Bracket in Fourier Space as 

(2.5) (/,9) = ^|£^(kxl)a„+,, 

k,l 


The corresponding infinite dimensional Poisson structure matrix is Jk,i = (k x 
l)aJk+i- Then (12.4p is a non-canonical Hamiltonian system with corresponding 
Hamiltonian 


( 2 . 6 ) 


E 

kGZ2\{0} 


W+kW-k 


- E 

kGZ2\{0} 


W+k^ 

lkp ■ 


The Hamiltonian is obtained from the Kinetic energy 



u|pdx 


J D4'dx. 


2.2. Galerkin Truncation. We now truncate to a finite mode approximation and 
study the spectrum of the equilibrium corresponding to D*. We will present two 
approaches to this: a Galerkin-style finite element truncation, and a more sophis¬ 
ticated Poisson structure truncation by Zeitlin. 

First consider the Galerkin-style truncation (see m)- Define the domain for our 
truncated Fourier modes 

(2.7) V =[-N,N]‘^ 

Now set (Uk = 0, d)k = 0 for all k /gD. Then the differential equations (12.41) 
define a finite set of ODEs, but not a Poisson system. 


2.3. Zeitlin’s Truncation. An alternative truncation is that described by [23] 
(see also Eq], M)- Restrict to the set of Fourier modes to 

(2.8) (4;k, k G D, 


and whenever a mode is referenced that is outside the domain 27 it is mapped back 
into 27. For this we have the notation k, which for any k denotes a mode k G 27 
for which the difference k — k = {2N -|- l)(a, 6)‘ for some integers a, b. 

Zeitlin gave the following Poisson bracket on the domain 27: 


(2.9) 


( 2 . 10 ) 


{f,9} = E 

k,iex) 


sin(ek x 1) df dg 


dujk dull *^+’’ 


Tk.i = - sin(ek x l)o 


where k, 1 G 27, and e = 
Hamiltonian 

( 2 . 11 ) 


e 

2-77 

2iV+l ■ 


k+1 


The corresponding truncation of (12)61) is the 




a2+ka2_k 


kGX>\{0} 
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where only the domain of summation has changed. 

The vector field under the Zeitlin truncation is thus given by 


( 2 . 12 ) 

(2.13) 


Wk 


(JV?^)k = ^ Jk.iVi/i 

IGD 

= i^sin(ekxl)caj^^jj^. 

^ ic n II 


The primary theoretical advantage of the Zeitlin truncation over the Galerkin 
truncation is that 27V +1 of the Casimirs present in the full system are preserved in 
the Zeitlin truncated system. The disadvantage of the Zeitlin truncation is that we 
must take some care when making arguments based on it (see Section |4]). The two 
truncations will be compared both numerically and analytically in later sections. 
For details of the construction and a description of the Casimirs see [Ml ES] . 


2.4. The Linearised system. For a fixed p G 2? \ {0} the equilibrium point of 
the PDF n* = 2F cos(p • x) is an equilibrium point of the truncated ODE given by 


(2.14) 


r F if 1 = ±p 
( 0 otherwise. 


As the original problem has symmetries x -O- —x, y -O- —y, and x -O- y, let 
P = {Pi:P 2 )^ with Pi > P 2 > 0 and pi > 0. 

The Jacobian of the Galerkin truncated vector field (12.4|) is 


(2.15) 


Ja,}a — 


0 ifa = borb = 0 



|b-ap 


b X a Wa-b- 


The Jacobian of the Zeitlin truncated vector field (12.131) is 


0 ifa = borb = 0 

(2-16) 


(a-b)' 


Evaluating these at the equilibrium ()2.14p gives the linearised systems 


(2.17) Ok = F 


1 


1 


|p|2 |k + p|2 

for the Galerkin truncation and 
F ( 1 


k X p Wk+p - F 


IpP Ik-P( 


1 


k X p cjk- 


(2.18) 


tuk = — 1 _ sin(e(k - p) x k) + pr sin(ep x k) 

"Vlk-PP IpI" 


k—p 


( _L.. sin(e(k + p) x k) + sin(ek x p) | w. 

2£V|k + p|2 IpP ' 


k+p 


for the Zeitlin truncation. 


2.5. Decoupling into Classes. The key observation is that iFk depends only 
Wk±p- Thus the linearised systems can be block-diagonalised. This block-diagonalisation 
is analogous to the construction in m- Following Li we call the individual blocks 
classes, which leads to the following definition; 
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Figure 2.1. The differential equations governing the set of 
Fourier Coefficients decouple into ‘classes’ when linearised. For 
a G Cjg, depends only on Wa+p and Wa-p- Extending this idea 
we get a subset of coefficients that only depend on each other, the 
class led by a. These coefficients all lie on a straight line with di¬ 
rection p. Classes that do not have an intersection with the disc 
indicated are stable. Clases that intersect the shaded region indi¬ 
cated at a lattice point have a pair of real eigenvalues (Theorem 
EH). Classes that intersect the disc but do not have intersect the 
shaded region at a lattice point lead to either a complex quadruplet 
or two pairs of real eigenvalues. 


Definition 2.1 (Classes). For some a G (and p fixed by the choice of equilib¬ 
rium), the class Ea C 2? is defined for the Galerkin truncation by 

(2.19) Ea = {a-hfcp G I fc G Z}. 
or equivalently for the Zeitlin truncation 

(2.20) E(^ = {a'T^ G P I fc G Z}. 

Figure [2T] illustrates this idea. Note that the Zeitlin truncated classes E(j ‘wrap 
around’ the domain D. 

m makes the analogous definition for the non-truncated system. In that pa¬ 
per, the classes are infinitely large, and there are infinitely many classes. In our 
definition, there are finitely many classes of finite size, which depend on the trun¬ 
cation size N. As V is finite, Ea and E(j are both finite. Write p = (pi,P2)^, and 
K = gcd(pi,p 2 )- Then 


2N + 1 


( 2 . 21 ) 


|Sa|< 


max(pi,p2) 
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and 

( 2 . 22 ) 


I Si 


2N + 1 

gcd(2A^ + 1, k) ' 


Note I Si I does not depend on a, and is odd for all choices of N and p. We use this 
fact many times later. For gcd(pi,p 2 ) = 1, |Si| = 2N + 1. 

For Ea we can make a “canonical” choice of a by selecting a in a length |p| 

II 

strip. For E^ a canonical choice for a is found by restricting a to be in a |p| by 
rectangle centred around 0 oriented so the sides of length |p| are parallel to p. For 
the Zeitlin truncation a unique choice is not possible if r > 1. See Section |4] and 
Id.ll for details. 

Fixing p and a canonical choice of a we now restrict our attention to the asso¬ 
ciated subsystems Ea and E^. Introduce new notation 


(2.23) 


UJk 


^a+kp if a -|- kpG'D, 
0 otherwise. 




^a+/cp’ 


(2.24) 



, — for k =G Z, 
|a-|- fcpP 


(2.25) 


a = Fa X p G R, a' 


F sin(£(a -t kp) x p) ^ 

e 


The value of the coefficient pk is related to the distance of the lattice point a -|- A:p 
from the boundary of the disc of radius |p|, where negative values of pk correspond 
to lattice points inside the disc. Note that a does not depend on k. Also note that 


(2.26) lim a' = lim a' = Fa x p = a. 

N^co £->-0 

Noting that sin(e(k -|- p) x k) = sin(ep x k) and rewriting (I2.17|) and (12.181) with 
the new notation gives a compact form of the linear systems 


(2.27) Wfe = a{pk+iu}k+i — Pfc-iWfe-i), 

(2-28) Wfe = a'{pk+iUJk+i - pk-iUJk_i). 

All indices in (12.281) are written modulo |E(j|. 

Let 1711,1712 G Z be such that a — mip G P, a — (mi -I- l)p and a-|-m 2 P G V, 


a + (m 2 -I- l)p /^'D. If we write w 

— (w—mi) 

mi + 1; 

1) Wo, W+1, 

UJrn2 

then cb = aAu3 where 







( ° 

+/3-mi + l 

0 

0 

0 

0 \ 


~P—mi 

0 

+P-mi+2 

0 

0 

0 


0 

~ P—mx + l 

0 

0 

0 

0 

(2.29) A = 

0 

0 

0 

0 

Tpm2 —1 

0 


0 

0 

0 

Pm2 — 2 

0 

“1“ Pm2 


V 0 

0 

0 

0 

Pm2 — 1 

0 / 
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If we write lo' = 

(^o> 

w'l, ..., 

l) 

T 

, then dj' 

= a'A'u}' 

where 


( ^ 

+Pi 

0 

0 

0 

~Pn-l^ 


-po 

0 

+P2 

0 ••• 

0 

0 


0 

-Pl 

0 

+P3 ■ ■ ■ 

0 

0 

(2.30) A' = 

0 

0 

-P2 

0 ••• 

0 

0 


0 

0 

0 

0 ••• 

0 

+Pn-1 


\+Pq 

0 

0 

0 ••• 

— Pn-2 

0 / 


If a = 0 or a is parallel to p, a = a' = 0. Thus the associated class only 
contributes zero eigenvalues and will not contribute to the linear instability of the 
system. We can thus ignore the classes with a = a' = Q. 

Note that A can be written a.s A = JS where 


(2.31) 


J = 


/ 0 +1 0 

-1 0 +1 

0-10 


0 

0 

0 


0\ 
0 
0 


fp 


— mi 



0 0 

P—mi + l 0 

0 P—mi+2 


0 \ 

0 

0 


0 0 0 

Vo 0 0 


0 +1 

-1 oj 


V 0 


0 0 


Pm2/ 


A similar construction can be made for A' = J' S'. As J and J' are skew- 
symmetric and S and S' are symmetric, these are (non-canonical) Hamiltonian 
systems. For both systems, % = From this it follows that if A is an 

eigenvalue of A or A' then — A, A and — A are also eigenvalues. Note that det( J') = 0 
as J' has odd size, and therefore J' is not symplectic with a one-dimensional kernel. 
J is symplectic if and only if jHal is even. 

We now focus on the behaviour of the eigenvalues of A and A' as a function of 
the pk values. Note that there is a symmetry in a. For every class Ea or E(j, the 
class E_a or generates the same set of eigenvalues. It is worth noting that 
a(—a) = —a(a) and a'(—a) = —a' (a), but as all eigenvalues occur in ± pairs, this 
does not affect the spectrum. Thus for the full system rather than a particular class 
all eigenvalues occur with even multiplicity. 


Definition 2.2 (The Unstable Disc). Introduce the disc Dp 
(2.32) Dp = {xeD| |x| < |p|}. 

This disc is shown in figure 12.11 A simple but important observation is 


Lemma 2.3. A lattice point is inside the unstable disc if and only if the corre¬ 
sponding p is negative: 

(2.33) a-l-fcpGDp Pk<0- 

Proof. This is true as < 0 if and only if |a -|- fcp| < |p| from (12.241) . which is 
exactly the condition that a -|- fcp G Dp. □ 


^As J' is circulant, one could write down its eigensystem explicitly by applying a discrete 
Fourier transform (see m) and therefore find a set of canonical coordinates (see m)- This will 
not be used in this paper, but may be of interest. 
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This is illustrated in figure [2T] The point a inside the disc corresponds to po < 0 
and the other points correspond to pk > 0. Also note that 

a + kp € dDp Pk = 0. 


3. Stability and Instability of Classes 

3.1. Stable Classes. The matrices A and A' are similar to skew-symmetric ma¬ 
trices by conjugation. 


(3.1) 


(3.2) 

TAT-^ 



^ V P—mi 

0 

0 

••• 0 \ 


0 

P—mi + 1 

0 

0 

T = 

0 

0 

y/ P—mi+2 

0 


V 0 

0 

0 



/ 

0 

“1” yj P—mi P—mi + 1 

0 

••• 0\ 


yj P—m\P—7rL\ + l 

0 

+ y/ P—mi + lP—7711+2 

... 0 


0 

y/ P—mi + lP—mi+2 

0 

... 0 

V 

0 

0 

0 

••• 0/ 


A very similar construction exists for A'. 

If p/c > 0 for all k, this transformation is real and thus A and A' are similar 
to real skew-symmetric matrices. Thus all eigenvalues are purely imaginary and 
A and A' can be diagonalised and so the class is linearly stable. By (I2.33P this 
condition is true exactly if 


(3.3) Eani:'p = 0. 

This is the finite-dimensional analogue of Li’s Unstable Disc Theorem (Theorem 
III.l) in [n], though the method of proof used in that paper is naturally very 
different. A discussion of the details such as choice of p and N required for (|3.3I1 
to hold follows in Section 2] Because of this result, only classes with a G Dp can 
contribute linear instability. Also a = 0 implies a = a' = 0 and so this class cannot 
contribute linear instability. 


3.2. Unstable Classes. For classes with a G Dp, there are two primary possibili¬ 
ties to consider: 

i) There is exactly one intersection between the class and the disc (ie, Ea fl Dp = 
{a}). This can only occur when a is chosen to be in the shaded area indicated 
in figure [2A] 

ii) There are exactly two consecutive intersections between the class and the disc 
(ie, 

Ea n Dp = {a, a -I- p} or Ea H Dp = {a, a — p}). This occur when a G Dp is 
chosen outside the shaded area indicated in figure [2A] 

For the Zeitlin style truncation, there is also a third possibility we must consider: 

iii) There are at least two non-consecutive intersections between the class and 
the disc (ie, a G E^ fl Dp and a -f A:p G Ea D Dp for some k yf —1,0,1). 
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Note that points on the boundary are treated as being outside the disc. Also note 
that it is not possible for three consecutive lattice points in a class to be in the 
unstable disc. If a, a — p and a + p were are all in Dp, they would lie along a 
diameter as Dp has diameter 2|p| and the distance from a — p to a + p is 2|p|. 
Therefore a = (0, 0) and a ± p £ dDp. This is the only possibility to have three 
consecutive lattice point in Dp, and hence Dp can at most contain two consecutive 
lattice points. Figure I^TT] makes this idea clear. 

From our numerical and analytical results we can categorise the spectrum of the 
class in these three cases: 

i) The spectrum has a single pair of real eigenvalues and all other eigenvalues on 
the imaginary axis. This is proved in Theorem 13.II 

ii) The spectrum typically corresponds to a quartet of complex eigenvalues ±a ± 
(3i, and all other eigenvalues on the imaginary axis. It can also correspond to 
two pairs of real eigenvalues, though seems to be less common. 

iii) This corresponds to the class ‘wrapping around’ the truncated domain of lattice 
points and intersecting the disc again; see 14.11 The spectrum is a combination 
of case (i) and case (ii) according to how successive intersections with the disc 
occur. 


This last case is atypical and does not occur with the Galerkin truncation. Usually 
this case can be avoided by a proper choice of N, however when both entries of p 
are even, it cannot be avoided. This is discussed in detail in Section |4l particularly 
Lemma 14.11 

All our numerical evidence is consistent with the result in |16j that the number of 
eigenvalues with non-zero real part is < 2|I?p| (twice the number of interior lattice 
points in the unstable disc), and our observation is that this is the exact number of 
hyperbolic eigenvalues. 

For case (i). Theorem 13.11 proves that this case always leads to non-zero real 
eigenvalues. For some p and a leading to case (i). Theorem [3^ describes an explicit 
lower bound which is independent of the truncation size N for the real eigenvalues. 
This agrees with numerically observed results, eg the e = 0 inviscid result shown in 
Figure 2 of m- It should be noted that our methods do not preclude the possibility 
that there are other eigenvalues with non-zero real part unaccounted for; we simply 
assert that there is at least one eigenvalue with positive real part. This together 
with the results from [10] , [16] , and [22] is sufficient to conclude nonlinear instability 
for the whole system. 

For the following section, we consider a general set of parameters (uq, oi,..., a„-i) 
instead of pk- Introduce a tridiagonal matrix 


(3.4) 



^a+1 

0 

0 

0 

0 \ 

Uq, 

0 

aa+2 

0 

0 

0 

0 

— Oa+l 

0 


0 

0 

0 

0 

— Cla+2 

0 

0 

0 

0 

0 

0 

Uq,-|_3 

0 

0 

0 

0 

0 

0 

0 

0/3 

V 0 

0 

0 

0 

1 

1 

0/ 
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and its characteristic polynomial 

(3.5) Ti{x)=det{xI-TP) 

for some integers 0<a</3<n—1. Then T£{x) can be recursively defined by 
expansion from top left to bottom right 

Ta{x) = 1 , Ta^^ix) =x'^ + Qa+iaa, 

(3.6) r£(x) = xr£~^{x) + apap-iT^-^ix). 
or by expansion from bottom right to top left 

Tp{x) = l, Tp_i{x) = x^ + ap-iap, 

(3.7) Ti{x) = xT^^^{x) + a„a„+i7)f+2(2;)- 
Note that 


(3.8) Ofc > 0 for all a < fc < /3 => T^ix) > 0 for all a: > 0. 

This can be seen by the recursive definitions; all terms are positive. 0 
The following is also useful: 


(3.9) 


r£io) 


nLo 

0 


ak 


if /3 — a is odd, 
if /3 — a is even. 



d a 


fo 


if /3 — a is odd. 

(3.10) 


31 — 0 



if /3 — a is even. 


These can be proved by simple induction arguments. 
Introduce similar notation for (|2.30|) 


/ 0 

+Ol 

0 

0 ••• 

0 

1 ' 

-ao 

0 

+02 

0 

0 

0 

0 

-«! 

0 

+03 • • • 

0 

0 

0 

0 

-02 

0 

0 

0 

0 

0 

0 

0 ••• 

0 

+On-l 

V + dO 

0 

0 

0 ••• 

— an-2 

0 / 


with the characteristic polynomial 

(3.12) ^(a:) = det(x/ — A). 

Then, if n is odd, 

A{x) = xT^^-\x) + an-iar,-2V~\x) + a„_iaori"-"(x) 

(3.13) =7^-\x)+aQan-iT{^-^{x) 

This can be demonstrated by expanding by minors along the last row and column 
of xl — A. Recall that n is odd for the relevant problem (n = |Ea| from equation 

UM)- 


^Note that the 113.611 and 113.711 satisfy the condition for Favard’s theorem (see 0) • Thus the 
polynomials 'T{x\ ai, ...aj) are orthogonal for j = 1, 2, 3,... with respect to an inner product with 
some weight function (see [23]). However, as the a/c terms may be negative this weight function 
will not always be positive. This will not be used here, but may be useful in future work for 
describing the imaginary part of the spectrum. 
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We first show that in case (i) there is some non-zero real eigenvalue. Because 
of the Hamiltonian nature of the system this means there is a plus/minus pair of 
eigenvalues and so there is linear instability. This is then extended to show that 
under certain conditions there is pair of real eigenvalues with an explicit lower 
bound independent of N. 

Theorem 3.1 (Real Eigenvalues in case (i)). If po < 0, and pk > 0 for all k = 
1 , 2,..., n — 1 and pk = 0 for at most one of k = 1 , 2 , ..., n — 1 , then for sufficiently 
large N, (I2.30|) has a non-zero real eigenvalue. 

Proof. The characteristic polynomial A for odd n has leading term x" and constant 
term 0. By combining (13.131) and p.lOl) (noting that n is odd, so n — 1 and n — 3 
are even), the linear coefficient is given by 


(3.14) 

(3.15) 


dx 




x—0 


dx 


aottn-l ^(x) 


n—1 


0 

n-2 


E n ' 

fc=0 \j=0; 


E n ' 

\i=l; i^l+2fc 


n-1 / n-1 \ 

=E n 

fc=0 \i=0; i/fc j 





Note that (13.151) is only valid for Ok 0 for all k, but (13.141) is always valid 0 . Now 
let Ok = Pk- First assume pj, > 0 for all k As pk ^ > 0 as |a -|- A:p| —>• oo 

and the size n of the classes grows linearly with N, ^ > 0 for 

sufficiently large N (where pk = ak). However 0^=0 < 0 Eis oq <0, Oj > 0 for 

all j = 1 ,.., n — 1 , and hence the linear coefficient of the characteristic polynomial 
is less than zero. 

If Pfc = 0 for exactly one k, then (13.141) consists of only one term, OjE j/fe Pp 
This is less than zero as po < 0 and pj > 0 for all j ^ 0, fc. 

If Pfe = 0 for more than one value of k, then the linear term is zero and we cannot 
apply this argument. 

As the constant term is zero, and the linear term is non-zero, then the lowest 
order non-zero coefficient of the polynomial is negative. 

We now argue by contradiction. Assume all roots of the polynomial are imag¬ 
inary (say iiOk) or complex ( 7 ^ -\-i5j) or zero. Then because eigenvalues occur in 


^The expression in equation II3.14II is tiie so-called n-1®* elementary symmetric polynomial in 
the variables aj 
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positive and negative pairs as well as conjugate pairs the polynomial has the form 

(3.16) 

A{x) = 

712 713 

x"! (x - iuJk){x + iujk) [{x - - idj){x - -fj + i6j){x + -fj - i6j){x + 7 ^ + iSj)] 

k=l j=l 

(3.17) 

712 713 

= n(x" + ujD {x* - 2x^(n] - (5|) + ( 7 I + (5|)2). 

fc=i j=i 

The lowest order non-zero coefficient (of x^^) is YYjLi {{ij + ^f)^) > 0- 

Thus by contradiction there must be some real eigenvalue, which will occur in a 
plus and minus pair. 

□ 

For the Galerkin case (EH), although the result is the same we cannot apply the 
same proof. Since we were not able to show that these eigenvalues remain bounded 
away from the imaginary axis when N ^ 00 we now show in a different way that 
for both truncations there are eigenvalues whose real part is bounded away from 
zero for iV —^ 00 . These proofs have stricter conditions on the pk then Theorem l3.ll 
but we will see that they can be met for most p. 

Lemma 3.2 (Lower Bound for Real Eigenvalue (Zeitlin)). If gq < 0, and Ofc > 0 
for all fc yf 0 , and qq + a 2 < 0 , then 

(3.18) .4(\/—ai(ao + 02 )) < 0. 

Proof. By expanding (13.131) using p.7p . 

(3.19) A{x) ={x^ + (aoai -I- aia 2 )x)Tff~^{x) 

+ {x^ + aoai)a 2 afTi~^ {x) 

+ aoan—iTi ix). 

Thus 

(3.20) A{\J—ai{ao + 02 )) = — aia\afTff ^(v^—ai(ao + 02 )) 

+ ooa„_i77* ^(-v/— 0 . 2(01 -I- oa))- 

As Oi, 02 , 03 ,Oji-i > 0 , Oq < 0 and the T terms are positive by (13. 8p the result 
follows. □ 

We make a similar case for the Galerkin truncation 

Lemma 3.3 (Lower Bound for Real Eigenvalue (Galerkin)). If fi > 0, gq < 0, and 
Ofc > 0 for all k 0, and ao + 02 < 0, then 

(3.21) < 0. 

Proof. Begin by noting that 

(3.22) 77(\/-oi(oo + a2)) > 0 
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for all 7 < 0 by (|3.8I1 . By expanding (13.7p . 

(3.23) Ta{x) = {x^ +aiaQ)T~^{x) + aoa-ixT~‘^{x) 

(3.24) — 0 'i{aQ + 02 )) = — 01027 )), ^(\/—ai(oo + 02 )) 

+ aoo_iv^-ai(ao + 02)7)7^('\/-oi(ao + 02 )). 

As Oo < 0, Ofc > 0 for all A: ^ 0, and T~^ ,T~^ take positive values for positive 
arguments, Ta{\/—ai{ao + 02 )) < 0 . 

Now similar to the proof of l3.21 

(3.25) Ta{\/-ai{ao + 02 )) =(x‘^ + 01(00 + a2))xT~^{\/-ai{ao + 02 )) 

+ (a:^ + a 2 ai)aoa-iT~'^{^-ai{ao + 02 )) 

= -aiOoa_i7)7^(-\/-ai(ao + 02 )). 

As oi, o_i > 0 , oo < 0 and the T terms are positive by (13.81) it follows that 

-oi(oo + 02 )) < 0 . 

Now if 7 > 2, we can make a recursive argument: 

(3.26) rj(x) = xT2-\x) + o^a^_i7J-2(cc). 

Now X = y^—ai{ao + 02 ) > 0 and a-yU^-i >0 (as 7 > 2). By inductive reasoning, 
as Ta < 0, < 0, then 7)7 < 0 for all 7 > 2. Therefore T7(-\/— Oi(oo + 02 )) < 

0. □ 

We can also make a similar construction if oq < 0, Ofc > 0 for all fc ^ 0 and 
Oo + o „_2 < 0. In this case, A(-\/—o„_i(ao + 0 ^- 2 )) < 0. 

Lemma 3.4. Assuming the same conditions as Lemma \3.‘A and Lemma \3.S\. there 
exists some 

x*,X 2 > \/-ai{ao + 02 ) 
such that T^{x\) = 0, ^(x^) = 0. 

Proof. The leading order term of A(x) is always x" regardless of oo,..., a„_i, and 
A{x) is real for real X. Thus lima,_>oo ^(x) >0. But bv Lemma lHT^ Al^/— oi (an + ao)) 
0 , and the result follows by the intermediate value theorem. 

The same argument can be applied to 7)f(x), with the equivalent result. □ 

We now turn our attention back to the context of our problem. 

Theorem 3.5 (Lower Bound for Real Eigenvalues in case (i)). // a S Dp and 
a + fcp fzDp (equivalently a + A:p fzDp) for fc = 1 , 2 ,.., n — 1 and 

(3.27) = \/—Pi(po + P 2 ) 

is real, there exist Ai, A 2 > A^ such that Ai is an eigenvalue of (|2.29p and X 2 is an 
eigenvalue of (j2.30l) . 

Similarly if A^ = —pn-i{po + Pn- 2 ) is real there exist Ai, A 2 > A^ such that Ai 

is an eigenvalue of (12.2911 and X 2 is an eigenvalue of (12.3011 . 


Proof. This follows from the previous three lemmas, letting pi = Oi and making 
note of (12.331) . □ 
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Section 0] clarifies under what conditions is real, and Theorem 13.51 holds. For 
example, when p = (1,1)^ then a = (0, 1)^ leads to a case (i) class, so Theorem l3.ll 
holds, but the reality conditions in Theorem 03] are not satisfied. It is also possible 
that = 0, so the real eigenvalue could be zero, so we also require that A^ > 0. A 
sufficient condition on a for this to hold is given in Lemma 14.21 


3.3. Associated Eigenvector. We must also consider whether this eigenvalue is 
“valid” in the sense that the corresponding eigenfunction in the full system is in 
fact a function in L^(T^), the square integrable function on the torus. For this we 
need to show that the Fourier coefficients are in In this section we are going to 
work with the infinite dimensional system. The reason for this is that although the 
truncations are useful in the calculation of eigenvalues, it is simpler to analyse the 
corresponding eigenvectors in the full system. Thus our approach is to use the finite 
dimensional Zeitlin truncation to compute the eigenvalues, then take the limit as 
A —>■ oo, and then study the decay of the corresponding eigenvector. 

Consider the linearised matrix for the full system 


(3.28) 


0 

Po 

0 

0 

P-1 

0 

Pi 

0 

0 

-po 

0 

P2 

0 

0 

-Pl 

0 




/ 


This is the limit in some sense of the matrices (I2.29|) and (12.3011 which corresponds 
to the limit of the Hamiltonian systems A = JS and A' = J'S' with Hamiltonian 
'^i^) = T,kPk^k- 

Lemma 3.6. The infinite dimensional linearised system given by (|3.28l) has a 
positive real eigenvalue A under the same conditions as Theorem \S.5[ Furthermore, 
the associated eigenvector is in P'. 


Proof. According to Theorem 13.51 there exist positive real eigenvalues with some 
lower bound (independent of N) of l|2.30|) (and (12.291) 1 for any N and either choice 
of truncation. By taking A —> oo we can conclude there exists some positive real 
eigenvalue A of (13.281) . 

Now consider an eigenvector associated with this eigenvalue, v = {...,v-i,vo,vi,V 2 , 
For this to correspond to a real eigenfunction of the full problem (that is, for the 
Fourier series to converge), we need these Fourier coefficients to decay sufficiently 
fast: they need to be a sequence in 

The entries of the (infinite dimensional) eigenvector Vk of (13.281) corresponding 
to eigenvalues A satisfy the recursion relation 


(3.29) 


Aufe = pk+iVk+i — pk-iVk-i- 


Since all pk ^ 0 this can be rewritten as 

(3.30) Vk+i = -^^Vk + 

Pk+l Pk+1 

Consider the limiting behaviour as fc —>■ oo. then pk In this limit 

solutions to (13.3011 behave like solutions to 

(3.31) Vk+i = A|ppufe + Vk-i, 


...f. 




























16 INSTABILITY OF EQUILIBRIA FOR THE 2D EULER EQUATIONS ON THE TORUS 


see, e.g., [H]- 

This linear recurrence has the general solution 
(3.32) vk = CimJ + C 2 M 2 

where Ci,C 2 G R are constants and are solutions to — A|pp^ —1 = 0 . 

Thus /ri/r 2 = —1 and without loss of generality |/ii| < 1, |/i 2 | > 1 (note that we 
cannot have |/ri| = \fi 2 \ = 1 as A|pp 7 ^ 0 ). 

Now as V is an eigenvector associated with a real eigenvalue, the span of the 
eigenvector is an invariant subspace of the Hamiltonian system with Hamiltonian 
7i{oj) = ~ ^(^) = e'^V. As the Hamiltonian 

is an integral of the motion, 'H(v) = 'H(e^‘v). By taking the limit t —)> — 00 , 
H(v) = ■H(O) = 0. 

Therefore, 


(3.33) 

^(v) = ^Pkvl = 

k 

(3.34) 

= -POVI 


k^O 


Now, if (72 7 ^ 0, 

(3.35) Pkvl Pk{C 2 plf -)> 00 , 

recalling that pk —>■ and pk > 0 for all fc 7 ^ 0. But |po| < 1 and Vk is finite, 

so there is a contradiction. Thus C 2 = 0 and Vk = Cip\ in the asymptotic limit, 
where \pi\ < 1. This is exponential decay, which is sufficient for the Fourier series 
to converge. 

Similarly for k ^ — 00 , the limiting behaviour is governed by 

(3.36) Vk-i =-X\p\'^Vk + Vk-i- 

Again, this means Vk is asymptotic to CiPi + C 2 P 2 Imi| < IM 2 I > 1. By 
the same argument as above, we can conclude that (7i = 0 and so Vk = C 2 P 2 as 
k ^ 00 . Thus the Fourier coefficients decay exponentially on both sides with |A:|, 
and hence v is in □ 


4. INSTABILITY OF EQUILIBRIA 

To prove instability of an equilibrium for a given p we need to find at least one 
unstable class Ea/S^. A necessary condition for a lattice point a to lead to an 
unstable class is to be inside the unstable disc. More precisely we desire a lattice 
point a that leads to a class of case (i) as this is the simplest situation for us to 
deal with. Theorem [33] asserts that there is a real eigenvalue with an explicit lower 
bound under some certain conditions. The goal now is to determine for which p 
there exists a lattice point a such that the conditions of Theorem 13.51 are satisfied. 

There are a number of other considerations when we take Zeitlin’s truncation. 
Although this preserves the geometric structure and the Casimirs of the original 
problem, it introduces a problem that is not present in m- This has already been 
mentioned in the classification of classes; it is the appearance of case (iii), in which a 
class intersects the unstable disc at non-consecutive points. This may occur because 
now we have periodic boundary conditions not only in physical space, but also in 
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Figure 4.1. For fixed p and a, there are some concerns with the 
wrapping of the Zeitlin truncation and the way this affects the 
intersection Ea fi Dp. These situations are discussed in Lemma 

O 

Left: the shortest distance between a and some non-consecutive 
a + /cp is at least ■ As we are interested in the limit TV —>■ oo, 
this distance can be made arbitrarily large, so that > 2|p|. 

This ensures that a and a + fcp cannot both be in the disc Dp. 
Right: The situation where there exists k G M. such that a + fcp 
lies on the line segment between a and a+p. This causes problems 
for our values of p. For a + fcp to lie at a lattice point on the line 
segment, gcd(pi,p 2 ) > 1- If gcd(pi,p 2 ) is odd, we can avoid this 
situation by choosing N per Equation USD; if Pi and p 2 are both 
even this situation is unavoidable. 


Fourier space. There are two distinct problems caused by this, as illustrated in 
figure 14.11 

The lattice points of a class lie on parallel line segments with direction vector p in 
the domain of Fourier modes. In the Zeitlin truncation there is more than one such 
line segment in the domain. The first problem appears when the distance between 
these line segments is so small that more than one line segment intersects the 
unstable disc. This can be fixed by making N sufficiently large. The second problem 
occurs when non-consecutive lattice points lie on the same line segment intersecting 
the unstable disc. If gcd(pi,p 2 ) is not even this can be fixed by choosing N per 
611). Note that for our purposes, gcd(pi,0) = pi and gcd(pi,p 2 ) = gcd(|pi|, |p 2 |)- 

Lemma 4.1 (Correct choices of N for the Zeitlin truncation). For all p = {pi,P 2 )'^ 
such that K = gcd(pi,p 2 ) is not even, there exists a sequence of N which increases 
without bound such that for all choices of a any two non-consecutive lattice points 
in cannot both be in the unstable disc. 


Proof. Let 
(4.1) 


(2TV-hl)K-l 


2 
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for S IN. Thus 2N + 1 = (2iV+ 1)k. If k is not even, then such an iV is a positive 
integer and thus a valid grid size. Select as a lower bound N > so that 

N > . p is fixed and finite so this lower bound is always finite. We can thus 

find an infinite sequence of N that increases without bound by letting N increase 
without bound. 

If X G Sa then x = a + A:p for some /c G IN. Thus x lies on the line parallel to 
the vector p that passes through the point a + Ax(2A^+1) for some Ax G (note 
that this point will be outside the domain D). 

Similarly if y G Sa then it lies on the line parallel to p that passes through 
a+ Ay(2iV + 1) for some Ay G Then the distance between these two lines is 
(4.2) 

|((a + Ax(2iV + 1)) - (a + Ay(2jV + 1))) x p| {2N + l)|(Ax - Ay) x p| 


(Ax — Ay) X p G so d = 0 or d > . If d > this corresponds to x 

and y lying on different line segments, as in figure HIT] Thus the distance between 
two points on different line segments is at least > 2|p| for our choice of N. 

If d = 0, then x and y must lie on the same line segment. Thus y lies on the 
line parallel to p passing through x. 

Write p = Kq so q = ( 91 , 92 )^ where gcd((j'i,q 2 ) = 1- Then as x, y G 
y = X + fcq for some fc G Z. But x, y G Ea so y = x + jp + {2N + 1) A for some 

A G 

Thus fcq = jp + (2iV + 1)A. So 


(4.3) fcq = jnq + (2N + 1)kA = ^(jq + (27V + 1)A). 


Thus fcq is divisible by k, but the elements of q both be divided by k by definition, 
so K\k. Thus y = x + K/3q for some /3 G Z, and so y = x+/3p. Then if |y —x| < 2|p| 
(the necessary condition for x, y both in Dp) this implies /3 = Oor/3 = ±l, and 
the result follows. □ 


The outstanding issue with Zeitlin’s truncation is that there is no appropriate 
choice of N when gcd(pi,p 2 ) is even. If k > 1 and gcd(K, 2A^ + 1) = 1, then p 
will generate all multiples of q due to the wrapping operation. Thus classes that 
intersect the disc can return after leaving the disc and intersect the disc again, 
breaking the assumption of l3.2l that pk < 0 for only one value of k. This behaviour 
continues for all values of N with gcd(K, 2N + 1) = 1. If it is even, this is true for 
any iV; if it is odd, we select appropriate N to avoid this. 

It is important to note that this is not an error per se, but merely a failure 
of Lemma KT\ for our proof. The wrapping of the Zeitlin truncation associates 
modes in an artificial way, but still generates correct results. For instance if p = 
( 6 , 2)^ = 2(3,1)^, the class led by a = (1, 6 )^ intersects the unstable disc again at 
(—2, 5)^ = (1, 6)^ — (3, 1)^ for any finite truncation size. If we compare the non- 
imaginary eigenvalues of the class led by a = (1,6)^ with those of the Galerkin- 
truncated systems for a = (1,6)^ and a = (—2,5)^, the same eigenvalues are 
generated, with similar convergence to Figure [521 

Fortunately, this problem does not arise with the Galerkin truncation as the 
wrapping operation is omitted, and so the proof goes through without issue. For 
both truncations it is still necessary to establish whether for a given p there is an 
a such that the conditions of Theorem 13.51 are met. We address that now. 
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Figure 4.2. If a is inside the dashed blue circle centred at the 
origin, then po + P2 < 0 and po + Pn -2 < 0. This circle has centre 
(0,0) and radius (-v/S — l)|p|. The shaded region here shows the 
overlap of this condition and the shaded region in figure 12.11 To 
show that there is at least one lattice point in the shaded regions 
we show the disc inscribed in this region (indicated by the small 
shaded red circles) has radius larger than i/\/2. 


Lemma 4.2. For all p = {pi,p 2 )^ except (1, 0 )^, ( 1 , 1 )^, ( 1 , 2 )^ (and permutations 
and sign changes thereof) there exists a choice of a. such that the reality conditions of 
Theorem lS. 51 are satisfied for an appropriate choice of N in the Galerkin truncation. 
Furthermore, if k = gcd(pi,p 2 ) is not even there is also a choice of a such that the 
conditions of Theorem 1,9.51 are satisfied for the Zeitlin truncation. 

Proof. For the bounds given in Theorems 13.31 and 13.21 to be real, positive, and 
hence a valid bound, we require Po < 0, pfc > 0 for all fc ^ 0, and Po + P 2 < 0 (or 
po +p „-2 < 0). 

If Po < 0, pi,p-i > 0, then |a| < |p| and |a± p| > |p|. This is true if and only 
if a is in the shaded region figure [2T] In the Galerkin truncation, this is sufficient 
to show that pk > 0 for all k 0. By Lemma 14.11 in the Zeitlin truncation we 
can find an unbounded sequence of choices of N such that this is sufficient to show 
fc ^ 0, pfe > 0. Thus for an appropriate choice of N we only need to prove that 
there exists an a such that |a± p| > |p|. 

As At = y —Pi(po + P 2 ) (or equivalently = y^—p-i^po + P- 2 )) is required to 
be real and non-zero, and p_i > 0 , then po + P 2 < 0 (equivalently po -|- p -2 < 0 ). 
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If |a| < (\/3- l)|p|, then |a± 2p| < |a| + 2|p| < (-/S + l)|p|. So 


(4.4) 


PO + P±2 < 


IPI' 
= 0 . 


1(^3-1)pP |(y3 + l)pP 


We thus need to show there exists some lattice point a such that |a| < (y3-l)|p| 
and |a ± p| > |p|. These two conditions are illustrated in figure 112] in the shaded 
region. Note that this is sufficient but not necessary for Theorem 13.51 to hold. 

The idea now is to specify that a is in the disc inscribed by the shaded region 
in figure 14.21 which we call Dc- This disc is tangent to the circles with radii |p| 
and centres ±p and the circle centred at the origin with radius (v^— l)|p|. It is 

a simple geometric exercise to show that such a circle has centre ^ 
radius - l) |p|. 

If a S Dc, it is outside the circle with centre — p and radius |p|. Thus a + p is 
outside Dp and pi > 0. Similarly, p_i >0. As Dc is inside the disc with centre 

origin and radius |p| clearly po < 0 and by the above po + p ±2 < 0. If 

we take a Zeitlin truncation, choose appropriate IV such that pfe > 0 for all P ^ 0 
by Lemma 4. Then the conditions of Theorem 13.51 are satisfied and the resulting 
bound A'l’ is real and positive. 

All that remains is to show that there exists an integer lattice point a G Dc- 
Any disc with a radius greater than must contain some integer lattice point (as 
it wholly contains a square of side length I). Thus if 


(4.5) 


|P|> 


V3 

V2(2 - 73) 


4.57. 


then the Dc has radius greater than 72 and such a lattice point exists. Note that 
this is a sufficient but not necessary condition on p . 

Checking the small number of p values with |p| < 4.57 and k odd there are ap¬ 
propriate lattice points a for most such p. The following table shows an appropriate 
value for a for most such p, and “None” where no such a exists. 


p 

(4,1)^ 

(3,3)^ 

(3,2)^ 

(3,1)"’ 

(3,0)^ 

(2,2)^ 

a 

(1,-2)^ 

(1,-1)^ 

(1,-2)^ 

(1,-2)^ 

(0,2)^ 

(-1,2)^ 


For reflections/rotations of these values of p the corresponding reflection/rotation 
of a is an appropriate choice. 

Thus for all p = (pi,P 2 )^ such that k = gcd(pi,p 2 ) except except (1,0)^, (1,1)^, (1, 2)^ 
and reflections and rotations of these, there is a choice of a so that the conditions 
of Theorem 13.51 is satisfied when N = fgj- jy > ^ 


Now we can combine the results about eigenvalues and eigenvectors from the 
last section and the conditions on p when they are applicable in our main 

Theorem 4.3. The steady state If* = acos(p • x) -|- /3sin(p • x) is nonlinearly 
unstable for all p = (pi,P 2 )^ except p = (±1,0) , p = (O, ±l) and possibly 

p= (±1,±1)^, p= (±2,±1)^, p= (±1,±2)^. 
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Proof. By Lemma [4.21 for all p except those listed above there exists some a such 
that po < 0 and po + P 2 <0 (or po + pn -2 < 0) for an appropriate choice of 
N. Thus by Theorem 13.51 there exists a real positive eigenvalue A. Moreover, the 
eigenvalue is greater than \/—pi(po + P 2 ) (or —pn-i{po + Pn- 2 )) which is both 
positive and independent of the choice of truncation size N. The truncation size 
N can be increased without bound, by Lemma 14.11 Hence there is a hyperbolic 
eigenvalue in the limit N ^ 00 and the spectrum of the PDE is unstable. Now 
recall that any steady state H* = a cos(p • x) + 2/3 sin(p • x) can be rewritten as 

= 2r cos(p • x) and so the full result follows. 

By Lemma [3.61 the eigenvector associated with the eigenvalue A is in The 
classes led by a and —a have the same eigenvalue, and the the corresponding eigen¬ 
vectors can be combined to construct coefficients Wk of a real eigenfunction 
corresponding to A. Since the eigenvectors v are in the periodic function Qx 
is in L^. Together with the result in m, which shows that the spectral mapping 
theorem holds, establishes linear instability. To conclude nonlinear instability we 
refer to the work of [10] and [22] . In [10] it was shown that sufficient conditions for 
nonlinear instability are linear instability together with a ‘spectral gap’ condition. 
In [22] it was shown that the essential spectrum of the linearised Euler operator in 
the cases we are considering is iR. Because of the presence of a point of discrete 
spectrum bounded away from the imaginary axis, we have a spectral gap, and hence 
nonlinear instability. 

□ 

Note that this does not preclude the possibility that the values of p listed as 
exceptions do not also lead to a linearly unstable steady state H*. In fact, for 
p = (l, 1 ) , ( 2 , l) and reflections/rotations thereof numerical results find non- 
zero non-imaginary eigenvalues. For p = (l,l) there is one complex quadruplet 
of eigenvalues, ±0.24822±0.35172/ to five decimal places, calculated with N = 1500 
and r = 1. For p = (2,1) there are two real pairs and two complex quadruplets 
of eigenvalues. 

Theorem 14.31 together with the numerical results mentioned and the spectral 
mapping theorem shown in [8] indicate that the only linearly stable equilibrium 
of type inn is p = (±1,±0) ,(±0, ± 1 ) , because only these have the single 
lattice point a = (0, 0)”^ inside the unstable disc. This exceptional case leads only 
to zero eigenvalues. All lattice points outside the unstable disc (including those 
on the boundary) do not contribute to instability (see section 3.1) and hence this 
equilibrium is spectrally stable. By [14] it is then linearly stable, and by [2] and 
m it is also Lyapunov stable. 

5. Some Numerical Results 

5.1. The Unstable Spectrum. The Zeitlin class decomposition means that we 
now typically compute the eigenvalues of {2N + 1) matrices, each of size {2N + 
1) X (2A^-|-1). Without the class decomposition, the eigenvalues of one {2N + 1)^ x 
{2N + 1)'^ matrix need to be computed. So the class decomposition results in an 
extremely significant saving of computation time. For a Galerkin truncation, this 
computational saving is even more pronounced. However, this is at the expense of 
accuracy (see figure [52]) ■ 
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All Eigenvalues forp=(5,3) 
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Figure 5.1. All eigenvalues for the case p = (5,3)^, F = i. 
The Zeitlin truncation is used with N = 200. Note there are 200 
(or 100 plus-minus pairs) eigenvalues with non-zero real part, and 
there are 100 interior lattice points in Dp \ {0}. This confirms 
the result from m- Of the non-imaginary eigenvalues, 56 are real 
and 144 are complex. The number of interior points that satisfy 
Pq < 0, pi,p-i > 0 is 24. All these points correspond to real pairs 
(two sets of 24 = 48). The other real pairs come from interior 
points with po < 0 and pi < 0 or p_i < 0. This usually creates 
a complex quadruplet but in a few cases corresponds to two real 
pairs instead. Increasing N does not change the number of non- 
imaginary eigenvalues. 


Figure 15.11 shows all the eigenvalues associated with a fixed value of p and N. 
There are exactly twice the number of interior lattice points in Dp \ 0. This agrees 
with the result in m that the discrete spectrum of the corresponding operator has 
at most 2|Z)p| — 2 non-imaginary eigenvalues. Our numerical results indicate that 
this bound is likely to be sharp; for all choices of p tested there is equality. 

Figure 15.21 shows the values at which the calculated eigenvalues converge as a 
function of the size N of our truncation domain D. Compared to the Galerkin trun¬ 
cation the eigenvalue converges for much smaller values of N when using Zeitlin’s 
method. 

Figure inSl shows the correspondence between the location of values of a and the 
types of eigenvalues of the class This corresponds to the results of Section |31 
Compare the positioning of the Fourier modes with figure I^TTI 

5.2. The Stable Spectrum. Figure 15^ shows the density of the imaginary parts 
of the spectrum for a = (—4,7)^, p = (7,5)^. There are also non-imaginary 
eigenvalues but these are not shown on the figure. 

For any e > 0, we can choose sufficiently large N so that there exists some a 
such that < e for all b G S(,. So the imaginary spectrum of this class can 




INSTABILITY OF EQUILIBRIA FOR THE 2D EULER EQUATIONS ON THE TORUS 23 


- Zeitlin Truncation 

- Galerkin Truncation 


’ I ^ ' 1 11 1 1 /' 


rriTjiTT^^ 

1 1 11 " 

i!" 


20 40 60 80 100 

N 


0.019 
0.018 
0.017 
0.016 
0.015 
: 0.014 
0.013 
0.012 
0.011 
0.01 
0.009 




- Zeitlin Truncation 

- Galerkin Truncation 


20 30 40 50 

N 


Figure 5.2. Numerically computed real eigenvalues vs Fourier 
mode domain size A^.For these figure a = (0,3)^ and p = (3,1)^. 
The red dashed lines shows the eigenvalues computed by the 
Galerkin truncation in equation p.29|l . and the black solid lines 
show the eigenvalues computed by the Zeitlin truncation in equa¬ 
tion (I2.30|) . For the left figure, the same truncation domain (12.7p is 
used for the Zeitlin truncation and the Galerkin truncation, mean¬ 
ing a larger matrix is computer for the Zeitlin truncation. For the 
right figure, a Galerkin truncation with 2N -(- 1 modes was chosen 
so that the same number of Fourier modes are included in both 
calculations. The convergence of the eigenvalue as a function of 
N computed with the Zeitlin truncation is significantly better in 
either case. These plots omit the factor of a/a' for clarity. 


be approximated by taking pk ~ The resulting matrix A from (12.301) is now 
circulant. A circulant matrix is diagonalised by a discrete Fourier transform (see 
m)- Thus the eigenvalues of A are then found to be 


(5.1) 


— 


2i 


■ sm 



for j = 0,..., n — 1 where j is the size of A. 


See [13] for details of this calculation. 

Thus the approximate imaginary spectrum of Ea for sufficiently large |a| lies in 
the interval ]^[~|ci|j |ci|] on the imaginary axis. Taking the limit N ^ oo (and so 
n —>■ oo), for each x € [0,1] there is a correspondence with an eigenvalue A^, where 
X = ^ sin“^ ( j. Differentiating this gives the density function 


(5.2) 


F{x) = 




That is, the proportion of the eigenvalues lying between cii and 02* on the 
imaginary axis is F{x)dx tor ci,C 2 G ■i^[—|a|, |a|]. This curve is also plotted in 
15.41 and it agrees well with the numerically calculated eigenvalues. This is surprising 
as the value of |a| is not particularly high. We can conclude that equation (15.211 
gives a reasonable approximation of the imaginary spectrum for many choices of a. 
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Eigenvalues with non-zero real part forp=(2,1) 
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Figure 5.3. Non-imaginary eigenvalues (left) and the correspond¬ 
ing lattice points a in the unstable disc (right) for the equilibrium 
with p = (2,1)^. At the top, the eigenvalues with non-zero real 
part for every unstable class Sa are shown. Eigenvalues with zero 
imaginary part are marked with a red x, complex eigenvalues are 
marked with a black -I-, and the zero eigenvalues are marked with 
a blue dot. In the bottom figure we see the values of a that corre¬ 
spond to these classes. The zero class led by a = (0,0)^ gives only 
zero eigenvalues. Compare the locations of the classes correspond¬ 
ing to real eigenvalues to the shaded region in figure [2Tj For these 
figures, F = |, so = cos(x • p), and a Zeitlin truncation is used 
for the approximation. 


[Hill] describe the essential spectrum of the linearised operator that coincides 
with our limit N ^ oo. The essential spectrum for the class led by a is given in 
that paper as 

(5.3) (less = 1^1], where ^ = |^(a x P)^. 

In the limit N ^ oo, axp, so —?> \f3\. Thus our approximation 

for large N reproduces the essential spectrum of a single class calculated in m- 
Note that this is the essential spectrum associated with a single subsystem of the 
linearised problem. The essential spectrum of the full system is the superposition 
of all these essential spectra. It was shown in that this is ill, which follows 
from considering the superposition of (15.31) for all possible values of a. 

6. Conclusion 

We have demonstrated the non-linear instability of the stationary solutions with 
vorticity ft* = 2r cos(p . x) in the Euler equations for almost all values of p such 
that p ^ (1,0)^ (or rotations and reflections thereof). We started out by using 
the Zeitlin truncation, and observe the numerical approximation of eigenvalues for 
finite N obtained from Zeitlin’s truncation converge much faster with N. However, 
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Figure 5.4. The density of the imaginary part of the spectrum, 
for the class p = (3,1), a = (1,-2), N = 1000, and T = 0.5. 
The bars show the normalized density of the imaginary eigenvalues 
computed for N = 1000. The thick red line shows the approximate 
density computed by taking the approximation pi~ —>■ . For this 
figure a Zeitlin truncation is used for the approximation. 


the Zeitlin truncation does not behave well when gcd(pi,p2) is even. Because of 
this, we developed much of our theory for both truncations. 

In addition, we have recreated and extended a number of results described by 
m, [IS], and [22]. Specifically, we have shown that the “unstable disc theorem” 
presented in m still holds true in the current context of finite dimensional approx¬ 
imation. Moreover, we have shown that for almost all p we can use the unstable 
disc to prove instability (as opposed to the existing stability results developed by 
Li and others). We have also numerically verified the bound on the number of 
non-imaginary eigenvalues and the essential spectrum of an individual class in |16] 
and [22]. We used very different approaches and arguments to those papers. 

There are obvious extensions of this work. The first is a complete description 
of the non-imaginary spectrum. This would require first showing that for two con¬ 
secutive negative values of p, the corresponding subsystem has four non-imaginary 
eigenvalues, either two real pairs or a complex quadruplet. This would be a step 
towards proving that the bound from [16] is sharp. Another extension would be 
to see if any of the methods used in this paper could be applied to more complex 
steady states, for instance fl* = sin(pia;) sin(p2y)- 

A similar analysis of the Euler fluid equations on a three dimensional torus would 
be interesting. A Zeitlin-style structure preserving truncation for the 3D case is not 
possible, as the Casimirs that make such a truncation useful are not present in the 
3D problemjj^. Similar stability results may still be possible using a Galerkin style 
truncation instead. 
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There is also a discussion of a structure-preserving truncation that includes a 
viscosity term in [26] . A comparison of Zeitlin’s truncation to standard truncations 
with a viscosity term included as in Figure [521 would be valuable. 

A significant extension of this material would be an application of the same 
methods to the Euler problem on a sphere. There similar structure preserving 
truncation also due to Zeitlin for the sphere [25] and so there is some hope of 
similar results in that setting. 
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